Opacity in the jet of 3C 309.1 



Eduardo Ros, Andrei P. Lobanov 

Max- Planck- Institut fiir Radioastronomie 
Contact author: Eduardo Ros, e-mail: ros@mpifr-bonn.mpg.de 

Abstract 

The "core" of a radio source is believed to mark the frequency-dependent location where the optical 
' depth to synchrotron self-absorption Ts « 1. The frequency dependence can be used for derive physical 

, conditions of the radio emitting region and the ambient environment near the central engine of the radio 

fSJ ' source. In order to test and improve the models to derive this information, we made multi-frequency 

dual-polarization observations of 3C 309.1 in 1998.6, phase-referenced to the QSOs S5 1448-1-76 and 
O '. 4C 72.20 (S5 1520+72) (4°26' and 1°49' away, respectively). We present here preliminary results from 

' these observations: total and polarized intensity maps, spectral information deduced from these images, 

and the relative position of 3C 309.1 with respect to S5 1448-1-76 at different frequencies. Finally, we 
discuss briefly the observed shift of the core position. 
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, 1. Introduction 

<N . . . . 

'sj- , Opacity in pc-scale jets The unresolved "core" of a compact extragalactic radio source is 

|0 I believed to mark the location where the optical depth to synchrotron self absorption Ts ~ 1. This 

position changes with observing frequency as i?core oc [1]. The power index kr depends on the 

shape of the electron energy spectrum and on the magnetic field and particle density distributions 

in the ultra-compact jet. Hence, by studying variations at kj. as a function of frequency, we may 

■ study the detailed physical conditions of the radio emitting region and the ambient environment of 
the source very near the central engine. Following [2] we can estimate basic physical parameters of 

O I the jet including luminosity, maximum brightness temperature, magnetic field in the jet, particle 
density, and the geometrical properties of the jet (i.e., the core location respect to the jet origin). 

■ To measure k,. a knowledge of the absolute position of the core (or at least the core offset 
between different frequencies) is needed. Hybrid maps in VLBI lack this positional information 
due to the use of closure-phase in the imaging process. The rigorous alignment of hybrid maps can 
be made by astrometric phase-referencing (e.g. [3, 4, 5, 6]). In sources with extended structure, 
an optically thin component can be used to align maps at different frequencies and then estimate 
the position of the core. 

The QSO 3C 309.1 The QSO 3C 309.1 (V=16.78, z=0.905^) is one of the most prominent 
compact steep spectrum (CSS) radio sources [7, 8]. Many of the CSS sources display polarized 
emission at cm- wavelengths. The ionized gas surrounding the jet disrupts its flow and is responsible 
of the complexity of the radio structures seen. There is evidence suggesting that 3C 309.1 is located 
at the center of a very massive cooling flow with M > 1000 Mq yr~^ within a radius of 11.5 h^^ kpc 

M- 

^This redshift corresponds to a linear scale of 5.60 pcmas"^ for _H"o=75 /i kms~^ Mpc~^ and qo—0.5. 
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VLBA observations of 3C 309.1 at 6 frequencies were used for determining the behavior of kr 
between 1.6 and 22 GHz [2]. To improve and extend this determination to higher frequencies, we 
observed 3C 309.1 using the VLBA at eight frequencies with dual polarization. In this contribution, 
we present a prehminary analysis of the new observations. 



2. Mapping analysis 

VLBA Observations We carried out VLBA multi-frequency, dual-polarization observations of 
3C 309.1 on July 19th and 23rd 1998 using all 10 VLBA antennas and observing at 1.5, 1.6, 2.3, 5, 
8.4, 15, 22, and 43 GHz. The data were correlated at the NRAO^ and processed using AIPS^ and 
DIFMAP [13]. DIFMAP was used for mapping the total intensity emission. The polarized intensity 
mapping and the phase-referencing analysis were carried out in AlPS. The 43 GHz data had to be 
discarded due to calibration and coherence problems. 



Total intensity mapping We applied the clean algorithm and self-calibration in DIFMAP to 
obtain the total intensity maps presented in Fig. 1 and described in Table 1. The resulting images 
show a structure very similar to the reported in earlier works (see [10] and references therein). The 
source shows a core-jet structure first oriented southward and turning later to the East at about 
60 mas from the core. 
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Figure 1. VLBA total intensity images of 3C 309.1. The synthesized interferometric beams are represented 
at the bottom left corner of each image. Contour levels are drawn at \/2 intervals. The left image follows 
the labelling convention from [12]. Map parameters (beam, total flux density, peak of brightness, lowest 
contour in map) are given in Table 1. 



Polarized intensity maps We applied the instrumental polarization calibration from the total 
intensity maps using the AlPS task LPCAL as described in [14]. We imaged the Stokes Q and U 

■^VLBA correlator, Array Operations Center, National Radio Astronomy Observatory (NRAO), Socorro, NM. 
^Astronomical Image Processing System, developed and maintained by the NRAO. 
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Table 1. Total intensity map parameters (Figs. 1 & 4) 
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a JVIinimum contour level in the figure. ^ Total flux density recovered in the map model. Corresponds to the B component. 
Corresponds to the A component. 



and produced images of the linearly polarized emission and the electric vector position angle. An 
exhaustive description of these results will be published elsewhere. We show an image of the linear 
polarization distribution at 1.5 GHz in Fig. 2. The core is unpolarized as in many QSOs. In the 
region South of B the electric vector is radial, suggesting a toroidal magnetic field viewed edgc-on. 
The degree of polarization is higher at the outer parts of the jet. 
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Figure 2. Polarized intensity image of 3C 309.1 at 1.505 GHz (A 21cm). Electric field vectors are shown. 



Spectral analysis The overall spectral index of 3C309.1 is —0.57 it 0.01. This result is obtained 
by adding together emission from all of the radio source structure which may have very different 
physical properties. 

To study spectral properties at different parts of the source, we mapped the radio source at 
all frequencies using natural weighting and very strong tapering (Gaussian function with half 
maximum at a distance of 33 MA). We convolved the CLEAN components with a circular beam 
of 4 mas in size, aligning the images on the peak-of-brightness of component A. We show these 
images in Fig. 3 together with the spectra of components A and B and the total spectrum of the 
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VLBI emission. The turnover frequency for A is around 8.4 GHz, and is below 1.4 GHz for B. A 
linear regression to the points for B provides an overall spectral index of —0.67 it 0.04. 
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Figure 3. Maps of 3C 309.1 at 4 mas resolution, obtained using the natural weighting and Gaussian tapering 
with the half maximum at 33 MA. The contour levels are drawn at \/3 intervals. The lowest level is of 
3 mJy/beam 



3. Phase-referencing analysis 

The calibrators We used two position calibrators for 3C 309.1: 

• 4C 72.20 is a QSO with z=0.799 and F=16.5, 1°49' East of the target source. It is a point-like 
source with an inverted spectrum. The imaging results are presented in Fig. 4. 

• S5 1448+76 is a flat spectrum, compact radio source with 2;=0.899 and y=20.0. It is 4°26' 
to the NW of 3C 309.1. It shows a faint jet to the NE at the lower frequencies, and it is also 
elongated in the East- West direction at the higher frequencies. The hybrid maps are shown 
in Fig. 4. 



The analysis We carried out the phase-referencing analysis in AiPS. We solved for the phase, 
delay and phase-rate for 3C 309.1, using the total intensity maps as input (dividing the (n, ■u)-data 
by the clean model) and thus removing the effect of the source structure. We then interpolated 
the values fitted using the task CLCAL for the fainter phase-reference calibrators, S5 1442-1-76 and 
4C 72.20 (S5 1520+76). 
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Figure 4. Hybrid maps of the position calibrator sources 4C 72.20 (top) and S5 1448+76 (bottom). The 
map parameters are described in Table 1. 



After editing the data, we mapped the radio sources using the MPS task IMAGR with the same 
parameters as were used for the hybrid imaging in DIFMAP. The phase-referenced maps are shown 
in Fig. 5 and the corresponding map parameters are given in Table 2. We measured the positions 
of the brightness peaks, whose offsets from the coordinate origin correspond to the offsets from 
the nominal position of 3C 309.1 relative to the calibrators. The relative positions deduced from 
this procedure are presented in Table 3. 




Figure 5. Phase-reference maps of 4C 72.20 and S5 1448+76 obtained using the phase, delay and phase-rate 
solutions for 3C 309.1 as described in the text. The contours are 49, 69 and 98% of the peak of brightness 
for each map. The values of the brightness peaks and their ratio (in percentage) with respect to their peaks 
in the hybrid maps (Fig. 4) are given in Table 2. 

The catalogue position of 4C 72.20 used at the correlator in error of +200 in a and -280 mas 
in (5. This translates into an estimated uncertainty of ~6 mas, in our preliminary position deter- 
mination at each frequency, making this fraction of the data unusable for our purposes. A proper 
analysis, correcting for the wrong position of 4C 72.20 will be published elsewhere. No ionosphere 
corrections have been applied in the data analysis. At frequencies lower than 8.4 GHz, the iono- 
spheric dispersion may severely bias our results. The tropospheric delay, especially the wet part, 
affects the phase for the highest frequencies, where the size of the water particles in the atmosphere 
is comparable to the wavelength. 

Notice that the ratio between the peaks of brightness of the phase-referenced maps and the 
hybrid maps (4**^ and 7**^ columns in Table 2) is the highest at the intermediate frequencies, where 
the compromise between the ionospheric and the tropospheric effects is found. Even when the a 
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Table 2. Phase-referenced map parameters for 4C 72.20 and S5 1448+76 (Fig. 5) 
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Tabic 3. Relative right ascensions and declinations in J2000.0 coordinates of 3C 309.1 with respect to 
4C 72.20, and S5 1448+76, obtained via AiPS phase-referencing. Biases have not been corrected. 
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No phase-referencing detection. 
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priori position of 4C 72.20 is in error, its ratios are similar to the ones in S5 1448+76, probably 
because the former is ~3 times closer to 3C 309.1 than the latter. 

The error budget in the positions (uncertainties in Table 3) includes the following error terms: 
a priori coordinates of the source, determination of peak-of-brightness in the maps, polar motion 
(estimated error of Imas), UTl-UTC (10~^s), station coordinates (5 cm), troposphere, ionosphere 
(oc and problems in the AiPS phase connection. This constitutes a conservative estimate of 

the uncertainty. We consider thus the phase-referencing results with S5 1448+76 at the highest 
frequencies as correct (central panel of Fig. 6). 

Offset 3C 309.1 - S5 1448+76 
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Figure 6. Left and central panel: Relative position of 3C 309.1 with respect to S5 1448+76, obtained via 
AIPS phase- referencing as described in the text. The data at the lowest frequencies are dominated by the 
ionospheric effect and are invalid for our purposes. Right panel: Comparison between the two different 
methods presented in the text to determine the core offset at 8.4, 15 and 22 GHz: the distance A-B with 
respect to the phase referencing results. The values at 22 GHz are set as zero reference. The A-B offset 
is biased due to beam effects in the L-shaped B component. The astrometric results are biased by the 
unmodcUed ionosphere and a simphstic model for the troposphere. 



The core position We assign the peak offset in a to S5 1448+76 and in 6 to 3C 309.1. The 
declination offsets in 3C 309.1 between contiguous frequencies are shown in Table 4. The relative 
offsets are also plotted in Fig. 6, where the value at 22 GHz has been set to be zero. 

Table 4. Core shift in declination for 3C 309.1 at the higher frequencies (see right panel in Fig. 6). 

Distance A-B AIPS phase- 
Frequencies in maps referencing 
15 -8.4 GHz -130+140 //as 150+430 fias 
22 -15 GHz 130+140 /xas 350+380 /ias 



An alternative way to measure the core offset is to assume that the B component in the maps 
from Section 2 is optically thin and its peak is at the same position for all frequencies. These 
values are also presented in Table 4, and in the right panel of Fig. 6. 

The trend in the dependence of the core position with the frequency is different in both methods. 
The A-B separation measurements is apparently frequency-independent: the difference at the 
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beams at different frequencies may bias this result, since the structure of the B component is 
L-shaped and in the dcchnation coordinate it is more extended to the South. Our prehminary 
positional results from the AlPS astrometry at 8.4, 15 and 22 GHz (but not at 5 GHz) suggest that 
the peak-of-brigthness of the maps shifts closer to the jet basis (core at infinite frequency) at higher 
frequencies, being this jet basis to the North of the A feature. This would be the expected opacity 
shift produced by the synchrotron self-absorption in the jet. Assuming that kr=l (self-absorbed 
core, [15]) at 8.4 GHz and that i?corc ^ u~^f^'' , the astrometric results provide values of k,- = l.l^f).b 
at 15 GHz and A:r=0.9±0.6 at 22 GHz. The big uncertainties do not permit to draw any conclusions 
about the physical parameters of the jet from kr at the present status of the analysis. A detailed 
analysis with the final, unbiased astrometric results will be published elsewhere. 

4. Summary 

We have presented preliminary results from a detailed multi-frequency study of the QSO 
3C 309.1 based on the VLBA observations made in mid 1998. We find a curved jet extending 
up to 100 mas to the East at low frequencies with two main components, A and B. The A com- 
ponent has a turnover frequency around 8.4 GHz and the B component is optically thin. The 
polarized intensity map at 1.5 GHz shows that the core is un polarized. In the region southern 
to B the electric field has a radial structure. The external parts of the jet have a high degree of 
polarization. A preliminary astrometric analysis provides a determination of the core position at 
diflcerent frequencies by phase-referencing to a nearby radio source QSO S5 1448+76. The changes 
at the core position with frequency suggest high opacity close to the core caused by synchrotron 
self-absorption. Dtic to the big uncertainties we cannot make any assert about the value of kr at 
high frequencies. An exhaustive analysis including ionospheric and tropospheric bias removal and 
physical modeling of the source will be presented in a forthcoming paper. 
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